gusucode.com > 《MATLAB图像与视频处理实用案例详解》代码 > 《MATLAB图像与视频处理实用案例详解》代码/第 16 章 基于Hu不变矩的图像检索技术/Compute_SimMat.m
function value = Compute_SimMat(I, J) [r, c] = size(I); A = []; for i = 1:r for j = 1:r M1 = (I(i, 2) * sin(I(i, 1)) - J(j, 2) * sin(J(j, 1)))^2; M2 = (I(i, 2) * cos(I(i, 1)) - J(j, 2) * cos(J(j, 1)))^2; M3 = (I(i, 3) - J(j, 3))^2; M0 = sqrt(M1 + M2 + M3); A(i, j) = 1 - (M0/sqrt(5)); end end value = A;